The symmetry of the superconducting order parameter in PuCoGas 
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The symmetry of the superconducting order parameter in single-crystalline PuCoGas (Tc = 
18.5 K) is investigated via zero- and transverse- field muon spin relaxation (/^SR) measurements, 
probing the possible existence of orbital and/or spin moments (time reversal-symmetry violation 
TRV) associated with the superconducting phase and the in-plane magnetic-field penetration depth 
A(r) in the mixed state, respectively. We find no evidence for TRV, and show that the superfiuid 
density, or alternatively, AA(r) = \{T) - A(0), are oc T for T/Tc < 0.5. Taken together these 
measurements are consistent with an even-parity (pseudo-spin singlet), d-wave pairing state. 



PACS numbers: 76.75.-(-i, 74.25.-q, 71.27.+a 

Because f-electrons may exist at the boundary between 
localization and itinerancy, f-electron compounds present 
an extreme challenge for many-body physics. For exam- 
ple, Pu metal sits midway between 'itinerant' Np and 
'localized' Am Q, while some theories for the 5-phase 
of Pu postulate that only one out of five of its 5-f elec- 
trons are delocalized |3| . Delocalization is a key ingredi- 
ent for heavy fermion behavior, which arises from strong 
f-electron conduction-electron hybridization. A new class 
of tetragonal heavy- fermion compounds Ce„MIn3„-|_2j 
with n — 1,2, has been under considerable investiga- 
tion for the wide variety of magnetic and unconventional 
superconducting behaviors which it exhibits j3] ■ The su- 
perconductivity in these compounds occurs below ~ 2 K, 
arising out of a normal state in close proximity to an 
antiferromagnetic quantum critical point and possess- 
ing a large Sommerfeld constant (7 = 0.7 J/mol-K^ in 
Celrln5). The recent discovery 4\ of superconductivity 
with an order-of-magnitude higher critical temperature 
in the same structural class of materials, but possess- 
ing Pu instead of Ce, is consequently both exciting and 
significant. PuCoGas becomes superconducting below 
Tc = 18.5 K from a normal state with a relatively mod- 
est 7 = 77 mJ/mol-K^ 4]. Photoemission suggests the 
f-electrons in PuCoGas are between localized and itiner- 
ant li]. 

The superconducting pairing interaction in heavy 
fermion materials is believed to arise from antiferromag- 
netic spin fluctuations, and a similar pairing mechanism 
has beenpostulated for the light-mass cuprate supercon- 
ductors This led to the speculation that PuCoGas 
might be a bridging material between heavy fermions 
and the cuprates ^ . Recent electronic structure calcula- 
tions for PuCoGa5 show 2-dimensional Fermi surfaces 
which could support spin-fiuctuation-mediated supercon- 
ductivity and, hence, d-wave pairing However, key 
measurements to confront these ideas are lacking; hence, 
tests of the symmetry of the superconducting order pa- 



rameter in PuCoGas are of supreme importance. 

This Letter reports muon spin relaxation experiments 
(/iSR) in single-crystal PuCoGa5 designed to explore two 
aspects of its superconductivity: (1) a possible violation 
of time-reversal symmetry (TRV), indicated by small, 
spontaneous magnetic fields below Tc, and (2) the tem- 
perature dependence and magnitude of its superconduct- 
ing penetration depth A. To date no direct spectroscopic 
probes of these quantities have been reported. In these 
experiments spin-polarized positive muons (/x"*", S = 1/2) 
are stopped in the sample and precess in the local mag- 
netic field B until they undergo a parity-violating weak 
decay /i"*" — > e+ + I'c + (lifetime Tp=2.2/zs). The 
time-dependence of the polarization, monitored via the 
positron decay anisotropy, yields information about the 
local field inhomogeneity and fluctuation spectrum 

Samples were grown from excess Ga flux forming 
flat plates with the c-axis normal to the crystal face. Two 
crystals measuring ~ 5x6mm^ in total area and ~l/2mm 
thick were encapsulated in a 70/xm thick Kapton coating, 
and then attached and sealed under He atmosphere in- 
side a titanium cell having a 50/im Ti-foil beam window 
and attached to a continuous-flow He cold-flnger cryo- 
stat. The dual encapsulation was undertaken to prevent 
possible radioactive contamination. 

Conventional time-differential /iSR experiments were 
performed at the M20 channel at TRIUMF, Vancouver, 
Ganada, only 3 weeks after the crystals were prepared. 
Thus, there was no measurable degradation of Tc due to 
radiation damage from Pu decay (half-life ~ 2.4 x 10^ 
y). Surface muons (momenta ~ 29 MeV/c) were im- 
planted into the sample with their spins perpendicular 
to their momenta, in the plane of the sample face. Ap- 
proximately 1/3 of the muons stopped in the sample, the 
remaining fraction in the Ti sample holder. A negligible 
fraction stopped in the Ti window or the Kapton coating. 
The background signal obtained from an empty Ti holder 
in either zero applied fleld (ZF) or 600 Oe applied trans- 
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FIG. 1: A representative zero field spin relaxation spectrum 
at T—4 K. The solid line is a fit to a dynamic Gaussian Kubo- 
Toyabe relaxation function and a background from muons 
stopped in the Ti cell. 



verse to the muon spin (TF) was well characterized by 
a Gaussian relaxation function G'Ti(i) = exp(— l/2a^;t-^) 
with (TTi ~ 0.03/is~^. For ZF experiments the residual 
field was reduced to < 10 mOe using trim coils. 

Figure ^ shows a representative ZF spin relaxation 
spectrum at T—AK, corrected for differences in the 
positron counters' efficiencies and solid angles. All of 
the ZF spectra - above and below Tc - were very simi- 
lar. Between T = 3 — 45 K the overall relaxation is well 
described by the sum of GTi{t) and a dynamic Gaussian 
Kubo-Toyabe spin relaxation function 8] G/fT(^, A, r), 
with static rate A due to a distribution of Co and Ga 
nuclear dipolar fields and dynamic rate r. The latter is 
produced by slow field dynamics and/or muon motion, 
as also seen in CeCoIus and Celrlns 9J. In PuCoGas 
these fluctuations are sufficient to suppress the long-time 
tail of Gkt (which approaches l/3AoGKT{t — 0) when 
r = oo), but only slightly affect the initial Gaussian re- 
laxation. The temperature dependencies of r"^, relaxing 
amplitude and rate A are shown in Fig. |21 No tempera- 
ture dependence is seen between T = 3 — 45 K. 

Transverse field /xSR locally probes the magnetic field 
distribution in the vortex state of superconductors [lo|. 
which is characterized by A and the coherence length 
estimated to be ~ 2.1 nm from the temperature de- 
pendence of the upper critical field Hc2 near Tc 01 • Field- 
cooled TF measurements were carried out in Hq ~ 600 
Oe field applied parallel to the c-axis. This field is at 
least twice the lower critical field Hci ~ 300 Oe reported 
previously Q and 4 ~ 5 x that implied by our measure- 
ments of A(0) {Hci oc ln(A/^o)/A^), as discussed below. 
Precession spectra were fit to the sum of two terms, cor- 
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FIG. 2: Temperature dependencies of inverse correlation time 
1/r, amplitude and field width A obtained from fitting the 
ZF spectra to a dynamic Gaussian Kubo-Toyabe relaxation 
function. 



responding to muons stopping in the sample and Ti cell, 
respectively: 

^oGz(t) = Acos{ujt + (t))P{t) 

+ ^TiCOs(wTii + 0)exp(-(TTii^/2). (1) 

The temperature dependence of a for a Gaussian 
P{t) = exp(— (T^t'^/2) is shown in Fig. One sees that 
a is much larger than (Txi ~ 0.03/iS~^ at all measured 
temperatures; thus, the two precession signals in Eq. ^ 
were easily separated. Fig.|2K also shows that a increases 
sharply below = 18.5 K due to the increasing field in- 
homogeneity caused by the superconducting flux lattice, 
and increases linearly within the statistical errors below 
about 12 K. 

A Gaussian time distribution implies a Gaussian field 
distribution, whereas in a single crystal the field distri- 
bution is expected to be asymmetric The quality 
of the Gaussian fits is illustrated in Fig. 01 for data 
taken at T = 4K. Here the G^{t) data (Eq. [T} from 
t ~ 3— 10/iS were fit separately, and this long-time Ti sig- 
nal was then subtracted from the total spectrum, leaving 
only the signal from the sample in the superconducting 
state. One sees that the Gaussian form for P{t) gives 
a satisfactory fit. Thus, the absence of the expected 
asymmetric field distribution is probably due to the large 
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FIG. 3: a) Temperature dependence of fiSK Gaussian TF 
rate <j{T). (See Eq. (1) and discussion in text.) The average 
baseline (Jn above Tc is indicated by the dashed line; the solid 
straight line shows cr ~ T below 12 K. The inset shows A(T) 
obtained from the linewidth crv(r) and Eq. 0. The solid 
line is a linear fit between 3 — 12 K. b) Rates from Gaussian 
((Tv) and exponential (Av - normalized to ctv at 3.3 K) forms 
for P{t) in Eq. (1) have the same temperature dependence 
below Tc. The curve is the result for the two-fluid model (s- 
wave), with Tc = 18.5 K and ~ 1.15/is~^, determined to 
fit the data for T > 12 K. The straight dotted line illustrates 
av(T),Av(T) (X T for T < 12 K. 



K = A(0)/^o ~ 100 value, together with the possibility 
that the flux lattice is distorted somewhat by radiation- 
induced pinning centers [l]| . 

We now discuss the data with regard to the super- 
conducting and magnetic properties of PuCoGas, be- 
ginning with the ZF results. First, the absence of any 
temperature dependence below 45 K in the parameters 
shown in Fig. |21is strong evidence that there is no static 
electronic magnetism. The shape of the observed short- 
time Gaussian relaxation shown in Fig. ^ is consistent 
with static nuclear-dipolar broadening. (If this relax- 
ation were instead due to damped, coherent precession 
from magnetic order, we estimate an ordered moment 
< 5 X 10~^iiB, which is <C 0.68/LtB obtained from the sus- 
ceptibility Assuming that muons occupy the same 
sites as in CeRhlus j we estimate a ZF nuclear-dipolar 
linewidth of about 0.29 fis~^ for both of the hypothetical 
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FIG. 4: Muon TF precession signal at T = 4 K. The data 
from t — 3 — lOfis has been subtracted to show the inhomo- 
geneous relaxation produced by the fiux lattice in the super- 
conducting state. The solid line is a fit for a Gaussian P{t) 
in Eq. (1). 



(1/2,1/2,1/2) and (0,1/2,0) sites . This is about 25% 
smaller than observed, indicating a different site occu- 
pancy from CeRhlns. Note, however, that none of our 
conclusions concerning the superconducting state depend 
at all on knowing the muon site(s). 

In TF measurements the /iSR rate from the vortex lat- 
tice (Tv measures the width of the rms field distribution 
{AB^y. Ctv = 7^(AB2)1/2 ^ where 7^ = 8.516 x 10* s'^ 
T"-'^ is the muon gyromagnetic ratio. The (AB^) is re- 
lated to A through the expression (assuming a hexagonal 
vortex lattice and Hq ^ Hc2) 

Ctv cx (AB2)1/2 ^ (0.00371)1/2$o/A2 oc 77s/m* (2) 

where $0 — hc/2e — 2.07 • 10^^ G-cm^ is the magnetic 
flux quantum, ?7s the superfluid density and m* the in- 
plane effective mass [l3| . 

The temperature dependence of CTv is obtained 
by subtracting the temperature-averaged normal-state 
linewidth CTn from the a shown in Fig. CTv = ct^ — ct^, 
where CT„ = 0.20 /xs~^. This assumes that the linewidths 
CTv and CTn arise from independent sources, which they 
do. To check the effect of the form of P{t) in Eq. (1) on 
the temperature dependence of the penetration depth, 
the time spectra were also fit using P{t) = exp(— Ai), 
and an analogous subtraction procedure was used (suit- 
able for exponentials), namely Av = A — An. As seen in 
Fig-Eb the exponential relaxation rate Av (normalized at 
the lowest temperature to the Gaussian rate CTv) yields 
essentially the same temperature dependence as does CTv, 
with a linear increase below about 12 K. 

The inset to Fig. shows the temperature depen- 
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dence of A(T), extracted from Eq. |21 We see that the 
characteristic length scale for A is '--^ 500 lattice con- 
stants {a = b = 4.232 A 01). Thus, if the muon makes 
~ 5 random hops in the 12/is duration of our experi- 
ment (l/r « 0.4/is^^) it moves only ~ -^5 lattice spac- 
ings, and cannot affect the measurement of A. We ob- 
tain an extrapolated zero-temperature penetration depth 
A(0) = 241 nm. To obtain the Ginsburg-Landau pene- 
tration depth one must apply mean-free-path corrections: 
Act = A(0)(l-|-£n/^t,r)^^^. where Ztr is the transport mean 
free path [l^- We have estimated Ztr ~ 10 nm using the 
values of the resistivity at T = Tc, Sommerfeld constant, 
Tc and S and the formulas found in Ref. 0|, so 
that Xgl ~ 265 nm. We note that our direct measure- 
ment of A(0) is about 2x that reported from critical field 
data '3| , and that the magnitude of A can depend on the 
fit function used Thus, our principal result is the 

linear temperature dependence of rjs (Eq. (0)), or alterna- 
tively, the low-temperature linacr behavior of A(r). As 
shown in Fig. (inset) AA(T) = A(r) - A(0) oc T for 
T/Tc < 0.5, within the statistical errors. 

We now discuss the implications of our measurements 
for the symmetry of the superconducting order param- 
eter in PuCoGas. The ZF ^SR linewidth A(r) shows 
no change in magnitude below Tc. Thus, there is no 
evidence for a TRV superconducting order parameter, 
the signature for which is an increased linewidth below 
Tc arising from spin or orbital moments [isj . as found 
in (U,Th)Bei3 Ea, Sr2Ru04 El, and PrOs4Sbi2 [13. 
This means that the superconducting order parameter 
is a linear combination of basis functions for tetragonal 
symmetry with real coefhcients. 

In a superconductor whose electrons are paired in an 
L = orbital angular momentum state the /iSR rate CTv 
or Av is relatively temperature-independent below about 
T/Tc = 0.5, reflecting exponentially-activated quasipar- 
ticle excitations over a superconduct ing gap which is non- 
zero over the entire Fermi surface This is clearly 
not observed in PuCoGa s, a s seen in Fig. where the 
two-fluid approximation |l3l | for s-wave superconductiv- 
ity is plotted as the solid curve. Instead, the relaxation 
rate continues to increase linearly at low temperatures, 
yielding a low-temperature T-linear behavior for rjs and 
AA(T << Tc). This behavior is quite generally associ- 
ated with a line of nodes in the gap function for strong 
spin-orbit coupling (as expected for the heavy element 
Pu))!^ and local electrodynamics [i^l; e.g., AA(T) cx T 
for T* < T <C Tc, where T* is a cross-over temperature 
T* vf/ttX{0), below which T^-dependence is calculated 
to occur due to non-local effects |2l|. We estimate the 
Fermi velocity Vf oc. CoTc « 3 X 10^ cm/s [H and obtain 
T* « 1.8 K. This is consistent with our results. A linear 
T-dependence has been found in a variety of copper-oxide 
high-temperature superconductors 01 > and is associated 
with d-wave {L = 2) pairing (even parity, pseudo-spin 
singlet). In many heavy fermion superconductors d-wave 



pairing has also been established [2j|. This particular 
gap symmetry is strongly enhanced by a spin-fluctuation 
pairing mechanism 6, 23], as opposed to electron-phonon 
pairing which favors s-wave pairing. Thus, our ZF and 
TF /iSR studies in PuCoGas show a gap symmetry which 
is similar to both the heavy fermion and copper-oxides 
superconductors [^l , though PuCoGas is neither partic- 
ularly heavy nor is it an oxide. This behavior may be 
tied to the 2-dimensional character of its Fermi surface 
which can enhance low-frequency spin fluctuations. 
Future /xSR experiments are planned to investigate the 
field dependence of the superconducting symmetry of this 
material. 
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